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ABSTRACT 


The resonance expansion method, developed earlier to study 
the propagation of sound in rigid rectanglular ducts, has been 
applied to the case of slightly soft ducts. Expressions for 
the generation and decay of various harmonics is obtained. The 
effect of wall admittance is seen through a dissipation function 
in the system of nonlinear differential equations, governing the 
generation of harmonics. As the wall admittance increases, the 
resonance is reduced. For a given wall admittance this phenomenon 
is stronger at higher input intensities. Both the first and second 
order solutions are obtained and the results are extended to the 
case of ducts having mean flow. 


I. 


INTRODUCTION 


This report deals with the propagation of higher order 
mode/ finite amplitude sound in slightly soft rectangular 
ducts., carrying mean flow. This report is a comparison of a 
report which deals with the rigid wall boundary condition (1) . 

Conventional methods of analysis have failed to solve, satis- 
factorily both these problems. ■ 

It was shown in the case of the rigid rectangular duct 
that the sound wave undergoes a nonlinear resonance. This is 
due to the fact that the eigenvalues of the higher order modes 
are integral multiples of the eigenvalue of the first mode. 

In reference 1 it was shown that a new method could be used 
to analyze such a problem. This method was called resonance- 
expansion method. It will be seen that the resonance pheno- 
menon persists for a soft duct as long as the admittance of 
the duct is not very large. Larger the amplitude of the 
sound waves , greater is the extent of the admittance to which 
this phenomenon will persist. 

Previous researchers (2,3) have attempted a solution of 
the soft wall problem by a direct multiple-scale expansion 
method. Unfortunately, the solution obtained by this method 
is not valid as the linear attenuation coefficient approaches 
to the order of the perturbation parameter e, (based on ini- 
tial amplitude) . This is due to the result that, in such a 
case, the direct method leads to secular terms in the second 
order solution. This alters the accuracy of the modal ampli- 
tudes in the first order solution. 

There exists an intrisic need to have a unified theory 
which explains the curious behavior in the entire domain of 
wall admittance. It would be seen that the solution described 
in this report approaches the one for the rigid wall case as 
the admittance tends to zero and to that from a straight for- 
ward solution as the admittance increases. There is also a 
pragmatic need for Such a solution as it promises a practical 
way to reduce resonance by means of a relatively small amoimt 


of sound absorbent. 

It should be noted that in view of the equivalent impe- 
dance concept of Cremer (4), (extended to higher order modes 
by Vaidya (5) ) , effects of viscosity and heat conduction could 
be Iximped on the wall of a duct as a small admittance. Thus 
this report could as well be viewed as one describing sound 
propagation in a rigid rectangular duct, carrying viscous, 
heat conducting fluid. 
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II. FORMULATION 

The problem formulated herein is idealized to the extent 
that an acoustic field is considered within an infinite duct, 
therefore only outgoing waves exist. (see Figure 1) 

The following set of non-dimensional variables are in- 
troduced, where { ) indicated the dimensional quantity: 

X = x/a, z = z/a, t = tc^/a, p = = ^/ec^a (1) 

where (j) is the nondimensional velocity potential. Therefore 
we have 

V = V4) (2) 

In terms of the non-dimensional quantities, the governing 
non-linear equations for the inviscid polytropic gas are (1) (e 


is amplitude of input wave). 

= eNj^(^) + e^N2((|)) (3) 

Ctj)) = + 2V({>- (V<j)^) (4) 

N 2 (cf>) = (y-1) (V(J))2 V^4>/2 + (V4>*V) (V(})) V2 ’ (5) 

Where 

( ) = 3^ ( )/3x^ + 9^ ( )/9 z2 (6) 

V( ) = [9( )/ X] I + [9( )/ z] k 

[V( )]^ = [9( )/3x3' + [3( )/3z]V (8) 

The relation between the pressure and velocity is (1) 

Yp = {(1“Y) + E^V(i))2/2] 4 - (9) 
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The formulation of the problem is completed once the 
boundary conditions are descirbed. Here we use a modified 
form of the semi-empirical form proposed by Zorumski aind 
Parrott (6) , (It's also used by Nayfeh and Tsai (3)) that is 

p-p = (R(u) - ~) u at X = 1 (10) 

° w at 

where u = is the velocity in the x direction. The form 
of this equation differs from that used in reference 3 in 
that the equation reduced to the familiar form 

P“Pq = (R + ix) u (10a) 

in the linear monochromatic case. 

In particular the following forms for R and X could be 
considered: 

R(u) = Rg + R 2 (u/Cjj) ^ (11) 

X(“> = Xo + X2 

It will be seen later on in this analysis that R 2 and X 2 pl^-y 
no part up to the second order solution. Therefore ordinary 
linear boundary conditions would be a(dequate. However, since 
the nonlinear boundary conditions are found in the current 
literature, we have chosen to retain them and to demonstrate 
their relative lack of significance. 

Substiduting equations (11) and (12) into (10) and non- 
dimensionalizing the resulting equation leads to 

p-y"* - «tV0o% ’ ^ ; 
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Expanding the right-hand side of equation (9) into Taylor's 
series, the polynomial relation between p and ((f) ,<}) ) can be 
found. Therefore equation (13) becomes 


at X = 1, [<f)^ + 



- Xo/t^ 


+ 0(e^)]<|), 


3t 


Where 



P c 
o o 



Xo 



(14) 


(14a) 


For the wall of duct with small admittance, let 0(n^)~e and 
0 (Xq)~ 1. Using this relation equation (14) can be arranged as 

at X = 1, (1 - Xq/w 3/3t) (f)^+E (ti-l/e) + O(e^) = 0 (15) 


The condition at x = O is 


X = 0 

(for symmetric case) 

(16) 

(j) = 0 

(for anti- symmetric case) 

(17) 


For the report, we consider the solution based on (16). The 
initial condition (or the input condition) at z = 0 is 

(f) = cos m X cos (jj t + 0(e) (18) 

Where m = miT (m = 1, 2, 3, . .. ,.) . The correction, 0(e) , is found 
such that the slightly soft wall condition is satisfied. (see the 
second order solution, section 3.3). 
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III PERTURBATION SOLUTIONS 

To determine an approximate solution for this system, we 
use the perturbation expansion. Let 


<P = <P^ + + O(e^) 


(19) 


where cj)^ and (p 2 are the functions of Zq, z^, ...... t and x. 


The z.'s are defined 
1 


z. = e z 
1 


( 20 ) 


Zq is a fast scale characterizing the wave length while z^ (i 
1, 2, 3,,.,,.) is a slow scale characterizing the amplitude 
modulation (xef . 2) . 


3,1 The First Order System 

Substituting equation (19) into the governing system, the 
first order system is (let be denoted as z) 


- <p 


o, tt 


= 0 


at z^ = O, = Re [cos mx cos o)t] (22) 


at X = 0, 
at X = 1, 


(p , = 0 
o,x 


b Jt. 


Therefore the solution of 4 is 

o 

N 


03 

(l) 


Ip . , 4 . = O 

o,xt 


4>q - .J, a . (z. ) cos jmxcos[je + a . ( z, ) ] 

J J J 


(23) 

(24) 


( 25 ) 


where 


= 03^ - m^, e = ax - 03t 
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Note that the form of equation (25) is of the resonent 
expansion type. This form could be generalized for various 
types of inputs. As a simple example, consider 

1 j = 1 


A. (O) = (26) 

J 0 3^1 

(0) = O j = 1,2,3....... (27) 


It would be assumed, at first that the order of the phases 
a^ is less than or equal to s. In Appendix B, this restriction 
is removed. It follows from equations (23) and (24) that 

m = in IT in = 1, 2, 3, ...... , (28) 

The behavior of A ^ ' s and a^ ' s can be determined by the con- 
dition that the second order system should be without secular 
terms. 

3.2 Solvability Condition for A ^ 

The governing equation for second order is 
N 

V^4)., - - S 2ja A*, cosjmx sin(j0+a.) 

X 3 3 

N 

- Z 2ja A.a^ cosjmx cos(j0+a.) = N '(^ ) (29) 

3=1 ^ ^ 3 I o 

In view of the assumption stated above, the phases shift 
can be neglected with the possible exception of its derivative, 
aV . Equation (29) could be rewritten as 


N 

V^(|), - ({)■, .. = E 2jaA! cos jmx sin j9 
X X , uu j=J_ ^ 


N 


+ Z 2ja a!A. cos jmx cos j0 + N, (<() ) (30) 

j = 1 3 3 X o 


where 


N N 

N. C<{)^) = Z E F^(x)j£2 A.A. sin (j+il) 
° j=l 5.=! ^ ^ 


N N _ 

+ Z Z (x) A.A sin (j“£)9 

j=l ^=1 ^ 3 


(31) 


Fj^(x) = [-J 5 (Y“ 1 )w -a w] cos jmx cos £mx + m^w sin jmx sin £mx 

(32) 


F 2 (x) = [-^ (Y“1) w^-a^<ja] cos jmx cos £mx - m^to sin jmx sin Jlmx 

(33) 

The boundary conditions are 


at z = O, = 0 

at X O, O 

at X = 1, X " X /w (j)^ = 

o l,xt 

Let the solution of (j)^^ be of the form 




(36) 


N' ~ N 

^ (x,z ) sin j9 4- Z Y. (x, z ) cos j9 (37) 

~ J X -;_i 3 X 


j=l 


j=l 


Substitution of equation (37) into (30) , (34) , (35) , and 

«Si • ' 

(36) yields the differential systems for 'F^ and 'F^ . For 'F^ , 
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'}’'I + H'. = 2 ja A*, cos jmx + 

3 3 3 


2 j A. »Ap [-Y- COS jmx + Y„ cos (j- 2 Jl)mx] + 

£,=1 3”^ X, X z 


[“Y^ cos jmx + Y2] (j/ 2 )^ A?^2 


E j£,(j + £) Aj_j_^A^ [“Y^ cos jmx + Y2 cos(j+ 2 £)mx] 


where 


= 


(Y+l) wV4 


^^2 = 


(o[m^ - (Y~l)a)^ -a^]/4 


^1 = 


n = 


j /2 - 1 if j /2 is an integer 


(j~l )/2 if j /2 isn't an integer 


N - j 


if j /2 is an integer 


O if j /2 isn't an integer 

The boundary conditions are 


at X = 0 > 4 ' . =0 

3 


at X = 1 , >F ! - j V 'P ! = 
j -'^o 

For 'i'j the equation is 


- (gj^/e)cos jm A ! 




2ja ajAj cos jmx 


( 38 ) 

( 39 ) 

(40) 

(41) 
'(42) 

(43) 

(44) 

(45) 

(46) 
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The boundary conditions are 


at X = O, 'i'l ^ Q 


at X = 1, + jx 'F'. = O 

Eliminating (1) between equations (45) and (48) yields 


at X = 1, 'i'! = 


-(rij^/e)jw cos jm A^/d + j^XQ^) 


It is shown in the appendix A that the system for '1'^ 
(equations (38), (44), and (49)) has a solution, if and 

only if 


A'. + f^A. - 


^1 :^i 


j 3 j 2a [ 2 Jl(j-il) A. A + 

Jl=l ^ 


^ (j/2)^ A^/2 ^ ^ ^3+^) ^j+£^£ ] = O, j = 1,2/3,... 


f = iil— *^o^o^o 
j ea ^2 ^ j2 -2 


ef. = 
J 


p c R 
o o o 


k^a"- m^TT^ R + j^X^ 
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where k is the dimensional wave number, (so that w = kc^) . 

From the dimensional form of f ^ , it could be recognized that 

dA 

contribution to 1 due to absorption is the same as obtianed 

in the linear theory by Doak and Vaidya (7). For other har- 
monics, equation (53) can be rewritten as 


ef . = 


ka 


/ k^a^ - m^TT^ 


Kt(j) 


where 


Kt(j) = r 


p c R 

OOP 

+ j^X ^ 

o '^o 


(54) 


(55) 


is the admittance ratio for the j harmonic. For the first 
harmonic, its definition coincides with the coventional defini- 
tion. Other harmonics have their admittance ratio computed 
after enlarging the duct reactance by j ^ . 

Equations (50) through (54) were obtained by assuming that 
the resistance of the linear was very large. Similar analysis 
could be carried out by assuming that the reactance was very 
large. Detailed analysis shows that the form of equations (50) 
and (52) is valid as long as the addmittance ratio is not an 
order of magnitude greater than e. The analysis is also valid 
above cut off of the exciting mode. 

As the admittance ratio gets smaller and smaller , the duct 
approaches the hard wall conditions. In that case the equation 
(50) reduces to the one obtained by Wang and Vaidya (1) . There- 
fore, there is a continuous transition from the hard to the 
slightly soft case. As the softness increases the generation of 
harmonics slows down. In that case, depending Upon the type of 
boundary conditions, simple multiple scale expansion or a straight 
forward expansion is adequate. 


1 


12 


3.3 The Solution of the Second Order System 

Having obtained A ^ ' s we could proceed with the second 
order solution. The condition that this solution be valid 
would lead to an equation in a^ with dA^/dz, replaced by 

equation (50) and from equations (35), (38) and (49). 


jY2 ^1 

^.(XrZ,) = -(a/m)f.A. x sin jwx + — y [ E A. .A. cos(j-2£)mx + 
1 1 5 3 4-nj2 36 


Z cos (j+2£)mx + ^ (A ?^2 “ ^j/2,1^^ ^ =1,2/3,... (56) 

4^/““ X 


where 


'^j/2,1 


1 

O 


j = 2 
j 7^ 2 


(57) 


Substituting the form of into the condition (48) gives 

at x = 1, a XpfjAj cos jm (58) 

Hence, for the system of 'i'j (equations 46,47, and 58) the con- 
dition of solvability is 


"j = 3 Xo 


(59) 


Therefore 


a. = j X f • 2- 

3 v^o 3 1 


(60) 


For small value of z (i.e. not far away from the inlet) the order 
of hj is e. This agrees with the assumption we made at Section 

3.2. The solution of is 

■ ■ 3 ... 


'FjCx,z^) = (a/m)jxQ f ^A^ x sin jmx , j = 1, 2, 3, . . . (61) 
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From equations (37)/ (56) and (61) the complete solution 
of the second order could be written as 


N jY2 

(j), = 2 { - (a/m) f . A. X sin jmx + — s- 
^ j=l J J 4m^ 


3l 

2 A. A cos(j-2£)rnx 

I •• n ^ A/ 

3=1 


^2 

2 

j=l 




cos(j+2£) mx + 


n/2 (Aj^2 




sin j0 + 


N 

2 {(a/m)j X sin jmx cos je} (62) 

j=l ° ^ ^ 


This solution is a generalized form of the one obtained for 
the case of the rigid walled duct. Equation (62) gives the 
correction, 0(e), of the input condition. That is 



(a/m)f^ X sin mx (x^ cos o)t - sin cot) (63) 
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IV* NONLINEAR PROPAGATION IN 
PRESENCE OF MEAN FLOW 


The governing equation in presence of mean flow was 
derived in reference 1. 


yV^(}) - - (j)^^ = e Ni (c{)) + e^N2(4)) (64) 

( 1 ) 

whAT-o II = 1 _ T " 1 , and N, ((Ji) becomes 

2 ^ 

N^((|)) = (y- 1) (<{>t - Mip^) V^tp + 2V(p» (V(p^ - (4) 


Therefore equation (21) should be changed to 


^J^jOfZZ ^ ^jo ^ ^jOfZZ ^^'^jo,zt ‘^jo^tt 


j “1/2,3, 


Assume again that solutions are the form in equation 
Substituting that into equation (65) yields 


0 

(65) 
(25) . 


a - — ^ r -Mo) ± [y (o)^ + M^m^ - ym^)]^Q (66), 

y-M 

Where the positive sign is to be used for waves traveling in 
mean flow direction and the negative sign is to be used for 
waves traveling against mean flow direction. y is the coef- 
ficient that takes into account compressibility. Note that 
in the conventional convective wave equation y is assumed to 
be unity. The soltuions of cj)^ and are of the same form as 
that in the previous section, except for the following modifi 
cations: 
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= = % (y+ 1) (a^+m^) (w-Ma) 


( 67 ) 


^2 “ ^2M ” %[(3“Y)m^ - (1+Y)a^] (w-Ma) 


( 68 ) 



Yi^/[a(ii-M^) + wM] 


( 69 ) 





( 70 ) 



are, once again, to be obtained by solving the equation 

The mean flow changes the transmission parameter a and 
These changes alter the propagation. The change depends 


upon the mode, the frequency of operation, magnitude and direc 
tion of the mean flow. 


V. NUMERICAL RESULTS 


Equation C50) was solved using N = 25 , Y^/ot = 15,6, to = 

3.685, a = 1.926, and four cases of wall impedance parameters 
shown in the table below; 


No. 

p C /R e 
•^o 0 0 

0 0 

• \ 

1 

0.0 

0.0 

0.0 

2 

1.0 

1.0 

0.75 X 10 ”^ 

'3 

10.0 

1.0 

7.5 X 10 “^ 

4 

1.0 

0.1 

1.48 X 10 “^ 

5 

50.0 

1.0 

3.75 X 10 “^ 


In this table the values of K^were obtained for j=l and 

-3 . . 

assuming e = 1.5 x 10 . It represents the admittance ratio 

for the principal mode. 

Figure 2 shows the effect of nonlinearity on propagation 
of first harmonic. The curves by dot line are obtained from 
linear theory with the boundary conditions given above. Which 
diminish only by the absorption of the wall. The nonlinear 
theory (solid line) shows more reduction in amplitude due to 
the transfer of energy into higher modes. 

Figure 3 shows for various impedances. As the admittance 
ratio goes Up, so does the attenuation. Figure 4 and 5 show 
the second and the third harmonic respectively. 

Figure 6 shows the ratio of A 2 to A^^ as wall admittance ratio 
changes. It should be noted that as the wall gets softer the 
generation of the higher harmonics is slowed down. For a very 
soft wall, the resonance would disappear and a straight forward 
multiple scale prediction would be adequate. 


VII. DISCUSSION AND CONCLUSIONS 


This report has outlined a method to obtain the sound field 
in a slightly soft rectangular duct. It was shown in an earlier 
paper that the rigid rectangular duct demonstrates a resonance 
phenomenon. The method of resonance expansion was proposed to 
solve that problem. In this report, the results are extended to 
the case when the wall admittance is small. 

It is seen that the phenomenon of resonance persists. How- 
ever, it gets weaker as the wall gets softer. The analysis 
shows that the criteria of softness is through parameters / 

R e and e C /X e. This implies that at high sound intensities 
o o o o 

even a moderately soft wall would act as a hard one. These para- 
meters show the classical antithesis between nonlinear convective 
terms and the dissipative terms. When the nonlinear terms pre- 
dominate over the dissipative terms, the resonance persists. 

It is also to be noticed that the same parameters which con- 
trol nonlinearity for the rigid duct control nonlinearity in a 
soft duct. Thus the parameter Y^/oi shows that the resonance is 
stronger near the cut-off frequency. Effect of mean flow is 
again qualitatively similar to that for the hard wall. 

It is interesting that for the principal mode the dissipation 
function is same as the one Which could be obtained by linear con- 
siderations. For higher harmonics the function is lpv;er than 
would be expected by the linear theory. 

The results also show that a small amount of absorption re- 
dices the resonance and hinders the process of higher harmonic 
generation. In a generalized form, this conclusion could lead to 
some practical applications. (see for example Reference 8) 
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APPENDIX A : COMPATABILITY CONDITIONS 


For the non-homogeneous linear system 


L(u) = r 


U^(u) =u* (b^) = Yi 


i = 1,2 


(A.l) 


where L is a second order linear differential operator and 
is some algebric operator indicating the boundary condition. 
The necessary and sufficient condition that this system 
may have a solution is that every solution, v, of the homo- 
geneous adjoint system (defined below) 

L(v) = 0 


Vjl(v) = V* (1) = 0, V2(v) = V' (0) =0 (A. 2) 

Satisfies the relation (9) which could be obtained from the 
Green's formula. 


vrdx = -YiV(bj^) + 


(A. 3 ) 


In second order systems if L = P — y + (2P'-P,) + {P'*-P' + P~) 

^ o 1 o 1 2 ' 

where the prime denote derivation with respect to x, in the 

particular case described in this appendix L =L, that is the 

equation is self-adjoint* In any case, any given equation can 

be converted into the self-adjoint form , rd /. dv , 

L = - G] . (see 


p. 215 of Ref, 9) The boundary conditions for the adjoint equa- 
tion are obtained from V. and K as illustrated below. 

^ - 2 
Consider the particular system defined as _ d ^ .2^2 
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In confirmity with^the rest of the paper let the argument u 
of this operator be ^ and v, the argument of the adjoint opera- 
tor be $ j . 

Therefore let 


L(J.) = -^ = (?.) + jV (J ) 


(A. 4) 


= Y. (0) = Yi = 


(A. 5) 


= Y2 


(A.6) 


Y 2 = -(n^^/e) jw A. cos jm / (1+j^)^^) 


(A. 7) 


r = 2jaA! cos jmx + E j£(j-«,) A. -A. [-Y 

J 2,=1 * 


cos (j-2£)mx] 


3, 2 


+ n[-Y, cos jmx + Y-] (j/2)-^A. 


mm 

- s jJl(j+£) A.. -A. 1-Y, COS jmx + Y- cos (j+2£)mx] (A. 8) 

jj,ssl x» J- A 

Therefore, the adjoint system corresponding to (A. 4), (A. 5), and 
(A.6) is 


L(^*-) = + j^m^$. = 0 


(A. 9) 


Vi($j) 


( 1 ) 


(A. 10) 


V2($.) = 


$j(0) 


(A. 11) 
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The adjoint system has a nontrival solution that is 


= cos jmx 


By the relation (A. 3), we have 


{A. 12) 


[2jaAj cos jmx + E j)l ( j-il) Aj_^A^j^ (-Y^^ cos jmx + Y 2 cos ( j-2fi,) mx) 


3 .2 


+ cos jmx + Y 2 ) (j/2) Aj ^2 


- E j£(j+^) (-Yj^ cos jmx + Y 2 cos ( j+2Jl)mx) ]cos jmx dx 

Xf“l 


.2 2 , 


cos jm( ~)jo)Aj cos jm/(l+j ) (A. 13) 


After integration, it yields 






(A. 14) 


I 
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APPENDIX B: GENERAL EQUATIONS OF A^ AND 


Form of the equation (50) when the phase angles a^ may 
or may not be small . 

Equation (60) provides meaningful insight about the 
order of magnitude of the phases a ^ . Rewriting it in the 
form 


(« -A ) «>aPo iXo 1 

a! ^ " ■ - — R — 5 — R- — 

^ a />2-2 -2-2 (R +j y ) 
Aj/ k a -m Ti ' o ^o ' 


(B.1) 


Order of magnitude analysis on this equation shows that 
, . jX 

wnen ^ would be of the order smaller than e. 

R D 

o 

j X 

When _o _ ^ would be of the order of e in the near field of 
R ' T 
o • 

the solution. In the farfield solution, a modified form of the 
first order solution is required. 

Working out the procedure outlined in Appendix A, when a^'s 
are not assumed to be arbitrarily small, equations of compatability 
become 

d(A. cos a. ) 

~ az ^ 

5,=1 

A._jAj^ oos(a.-j^+ aj^) + § (j/2)2 A^/2 cos (2a./2) 


"^^iJUj+iJ,) Aj^^A^ cos,taj^^-aj^)] = 0, j = 1,2,3,.. 


(B.2) 


I 

t 


23 



Note that equation (B.2) changes into equation (50) as a^ 0 

Equation (B.3) changes into equation (60) as a^ -)■ 0. These equa- 
tions will be discussed in a coming paper (10) which gives detail 
about the coup lex sound wave input. 
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Figure 1 ; Sketch of the semi-- infinite two dimensional 
rectangular soft duct. 

Figure 2 : predicted by linear theory (dot line) and 

nonlinear theory (solid line) , The conditions 
of curve 2 and 4 are given at Section V. 
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Ratio of A 2 to A^ for different B.C. The values 
within i-r) are (Pq^o/^o^ ' X^/R^) respectively. 
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